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Abstract 

Some fundamental Nucleon-Nucleon interactions and their applications to finite 
nuclei are reviewed. Results for the few-body systems and from Shell-Model calcu- 
lations are discussed and compared to point out the advantages and disadvantages 
of the different Nucleon-Nucleon interactions. 



1 Introduction 

A major task in our microscopic study of nuclei is to describe nuclear proper- 
ties with fundamental nucleon-nucieon interactions. First among these prop- 
erties is the experimental binding energy (of nuclear few-body or many-body 
systems). The prediction of this property remains a goal for any microscopic 
nuclear structure model; it remains a severe test of the underlying dynamics. 
The Shell-Model, in addition, also aims at predicting the spectroscopic details 
unraveled by the experiment. 

The first step towards a fundamental description of the nuclear few-body or 
many-body systems using the nucleon degrees of freedom is the establishment 
of the nucleon-nucieon (NN) interaction derived from the underlying dynam- 
ics; this can take various forms. From a fundamental point of view, one uses 
the derived realistic NN potentials to directly solve Faddeev equations for 
the eigenvalues of the nuclear few-body systems, or to determine a model- 
space dependent microscopic nuclear interaction for the nuclear many-body 
systems (to be used in the nuclear Shell-Model or the Bruckner Hartree-Fock 
approach). For nuclear many-body systems, starting from a fundamental stand 
point, one can take into account the strongly repulsive bare NN interaction 
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by first selecting an effective model space and then constructing an effective 
Hamiltonian (nuclear G reaction matrix) to describe the complicated many- 
body bound states. 

Effective Hamiltonians can also be defined based on simple parameterization 
of the NN interaction, with parameters fitted to the few body systems or 
directly for the large systems. This alternative avoids the use of the bare NN 
potentials. Effective Hamiltonians can also be defined to account for some 
characteristic properties (for instance, pairing correlation and/or symmetry 
issues 0) of a many-body system. 

At times, the effective Hamiltonians are the result of global fitting 0; this 
leads to the empirical Shell-Model approach. Let us say at the outset that this 
approach has achieved a very successful unified description for most of s — d 
shell stable nuclei 0. Yet, the Shell-Model has yet to produce a deep micro- 
scopic understanding. Besides, we know that the huge unfeasible Shell-Model 
spaces in the traditional large-scale Shell-Model calculations make it difficult 
to perform systematic calculations for nuclei beyond those in the f-p shell 
(except for some heavier nuclei near close shells). These nuclei often require 
huge scale Shell-Model calculations involving the best recent algorithms to 
deal with the large model spaces, while so many two-body matrix elements 
(195 independent 2-body matrix elements for fp-shell) cannot globally be fit- 
ted as easily as in sd-shell. Beyond the f-p shell, some dramatic truncation of 
the Shell-Model is required to even render the effective Shell-Model possible 
(for instance, broken-pair Shell-Model ||). Therefore, the microscopic Shell- 
Model study with NN interactions is not only tackling a fundamental nuclear 
many-body problem, but it may provide helpful information for the empirical 
approach. 

In this report, we will analyze the physics implied by different NN interactions. 
We will carry out simple Shell-Model calculations for modest valence-particles 
system (more specifically, 18 0, 18 F and 18 Ne ), where exact calculations can 
eliminate some uncertainties, starting from various NN interactions (i.e, using 
the NN interactions as input to obtain Bruckner's G-matrix as microscopic 
effective two-body interactions). We argue that different NN forces can give 
substantial discrepancies in nuclear low-lying spectroscopy. The aim of our 
calculations is to provide us with some lessons of and insights about the mi- 
croscopic study of the nuclear dynamics. 

This report is organized as follows: In sect 2, we will briefly review some mod- 
ern NN interactions and the physics they encompass. Sect. 3 will review the 
present status of the few-body and many-body calculations based on realis- 
tic NN forces. In sect. 4, attention will be paid to the nuclear many-body 
systems. The process to get effective interactions from NN forces is briefly 
mentioned. In sect. 5, we will first mention the recently developed Drexel Uni- 
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versity Shell Model (DUSM) algorithm and code, then describe some of the 
results and follow up by a brief discussion of the Shell-Model calculations. 
Finally, comments will be given in sect. 6. 



2 Nucleon-Nucleon Interaction 



The description of the nucleon-nucleon interaction remains one of the most 
fundamental themes in low-energy hadron physics. The resulting force has im- 
portant consequences for nuclear physics. With the advent of quantum chro- 
modynamics (QCD), one would naturally hope to derive the iViV force which 
is based on the quark-gluon dynamics. Typically, under the assumption of 
spontaneous breaking of dynamical symmetry and with 't Hooft's large N c ex- 
pansion approximation, and by formally integrating out the quark and gluon 
fields, one can obtain an effective chiral lagrangian of low-energy hadron dy- 
namics, namely, a "derivation" of skyrme-like lagrangian in the low-energy 
regime of QCD ||. Recently, a derivation from the QCD lagrangian to a 
hadron dynamics, lagrangian, which contains QCD string structured mesons 
and baryons, has been achieved. Based on the derived hadron dynamics la- 
grangian, the pseudoscalar coupling constant was theoretically estimated ||. 
Such "derivations" provide a justification for the effective meson model of 
nucleon-nucleon interaction. 



Due to the non-perturbative nature of low-energy QCD, the description of the 
NN force is carried out in practice directly via the nucleon-meson process, 
i.e., using the meson degree of freedom to describe the interaction between 
two free nucleons. Using this meson based description, the nucleon-nucleon 
interactions follow from three different approaches: second order perturbation 
theory (e.g., Yukawa potential) 0, dispersion theory (e.g. Paris potential) ||, 
and the field-theoretical meson-exchange model for the NN interaction (e.g., 
Bonn potential || and Nijmegen potential ||10|| ). There are also earlier and 
more phenomenological approaches to the nuclear potential; in particular we 
quote the hard-core Hamada- Johnston (HJ) potential the Reid soft-core 



(RSC) Potential [12 and the Argonne t>i4 ||13|| . We will come back to these 



later. We mention that chiral symmetry perturbation theory has recently been 
used to explore certain specific terms in the NN interactions (for example, 
the three-nucleon force) 



The NN force can be effectively divided into three interaction ranges: long- 



range one-pion-exchange part (r ~ 2 fin), intermediate-range attractive part 



(1 ~ 2fm) and the short-range repulsive part (r ~ lfm). This division follows 
naturally from the physics involved at these various distances. The long-range 
part presents fewer ambiguities since it is dominanted by the one-pion ex- 
change (OPE) effect. In fact, various realistic potential models have very sim- 
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Fig. 1. Average S wave AW potentials. 



ilar long range part; they differ however at shorter distances, and in particular 
in their intermediate-range part, which are supposedly due to the contribu- 
tion from heavier mesons (or pion resonance) or two-pion exchanges (TPE). 
For instance, the Paris potential determines the TPE from nN and mi inter- 
actions by using dispersion relations, while the Bonn potentials include both 
single-meson exchange (jr, uj, 5) and the TPE contributions (nNA, nAA, 
Np couplings) as general meson-field description. In Fig.l, an average AW 
potential for S wave is shown. The short-range repulsive terms have up to 
now been treated phenomenologically in all potential models. Since the mid- 
dle 60 's, the hard core is replaced by a soft-core. These forms are demanded 
to fit the experimental nucleon-nucleon phase shifts. A theoretical derivation 
of form factor functions for the potential near the origin would require further 
information on the subnucleonic structure and even explicit modelling of the 
quark-gluon degrees of freedom, which is beyond our current understanding 
of quark physics. Quantitatively, the contribution of OPE is about 95% of the 
total potential contribution for many nuclear properties, as is the case for the 
deuteron; in a sense the terminology "realistic potential" infers the realistic 
treatment of OPE and TPE. 

More specifically, in the Bonn potentials, all the coupling constants are de- 
termined by fitting np scattering data and the deuteron properties. It was 
claimed recently that the potential models which fit better the np data give 
poorer agreement for the pp data [fTo|] . From deuteron properties (i.e., the 
simplest np system), we already know that generally the np components (i.e., 
mainly T=0 states) of the nucleon-nucleon interaction are more attractive than 
the pp or nn (i.e., mainly T=0 states) components. Thus fitting only to np 
scattering data will create a two-nucleon potential which could overestimate 
the attractive components of the AW force with respect to a unified np, pp 
and nn many-body system. Recently, the Nijmegen group are attempting to 
provide a potential model by fitting both the pp and np scattering data. The 
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work to derive an effective interaction with this interaction is still in progress. 



3 Few- Body and Many- Body Methods 



There are many available modern potentials which are based on the analysis 
of Nit scattering data within some given energy regimes (for example, up to 
350 MeV, the threshold of % production ) and the deuteron properties. These 
potentials are the starting points for a microscopic study of the nuclear dynam- 
ics. To this end, the few-body system is an ideal testing-ground of the physics 
considered. Nowadays the methods used to obtain the nuclear few-body solu- 
tions (i.e., in configuration-space or momentum-space) provide very consistent 
and nearly exact solutions for the bound-state problem. Therefore, the results 
of the few-body problems will more or less free one from approximations and 
concentrate on the study of the nucleon dynamics. 

Solving the Faddeev-Noyes equations for bound states with a given nucleon- 
nucleon potential should provide us with the correct binding energy for such 
systems. Hence to study three-body bound systems, such as 3 H and 3 He, are 
particularly important since they are the smallest ones beyond the deuteron 
for one to quantitatively say something about their properties and the first 
set of nuclei for one to investigate the importance of hree body interactions. 
This seemingly simple goal has been pursued by many nuclear physicists for 
many years, and yet to date, there is still no definitive answer (for a recent 
review, see |I5|). For example, in Table 1 we have presented the benchmark 
results, which explicitly do not include three-body forces, of the predicted 
3 H binding energies based upon various modern realistic potentials and as 
functions of the number of three-body channels. Except for the Bonn potential, 
the calculated binding energies based on the various NN potential models are 
about 1 MeV underbound (the experimental binding energy is 8.536 MeV for 
3 H). The nuclear rms charge radii computed with realistic NN potentials are 
consistently larger than the experimental values. This discrepancy between 
experiment and theory is unnerving. 

One plausible explanation for the above problem is the presence of a three body 
force. Table 1 shows that a Hamiltonian of Reid soft-core two-body force plus 
the Tucson-Melbourne two-pion exchange three-nucleon force (TM3N) Jl6| can 
indeed provide additional binding for the trinucleon system. However, these 



results are highly sensitive to the 7rNN form factor and its cutoff mass A ||17|| . 
For instance, in a 34-channel calculation, the RSC+TM3N result for cutoff 
mass A=5.8 m, is 8.86 (MeV), while the binding energies are 7.46 (MeV) and 
11.16 (MeV) for A=4.1 and 7.1 respectively. Therefore, the current status is 
that while the correct binding energy can be reproduced with the adjustment 
of the three-nucleon force, it is done so in a phenomenological manner. On a 
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Table 1 

3 H binding energy (in MeV). 
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more fundamental level, three-nucleon interactions have recently been studied 



by the chiral invariant effective lagrangians |18| . This may be an encouraging 
direction to fundamentally understand and determine the effective nuclear 
three-body force. 

It has also been argued that the aforementioned discrepancy between the 
experimental and calculated binding energy is due to the missing of underlying 
baryon-meson degree of freedoms. The A-isobar and pion degrees of freedom 



should play some role in the trinucleon binding [19|. After all, the nuclear 



dynamics is just a low-energy limit of the baryon-meson dynamics. Sauer et al 



20| estimates that the presence of a single A-isobar provides about 0.3 (MeV) 
binding and that the three-nucleon force with an intermediate A-isobar gives 
0.9 (MeV) additional binding, whereas 0.6 (MeV) is lost due to the two-nucleon 
dispersive effect. Finally, a single A excitation produces about 0.6 (MeV) 
binding energy, which is not adequate to match the correct binding energy. 
While such effects can contribute to the binding energies of the respective 
three body systems, they are not adequate to fully resolve the discrepencies. 
Another factor which may muddy the water further is relativistic corrections, 
namely whether or not nuclear spectroscopy is insensitive to the off-shell part 
of the NN interaction. 



The transition from a nuclear few-body system to a nuclear many-body sys- 
tem requires new methodologies to obtain a solution since the traditional few- 
body techniques no longer apply due to obvious computational difficulties. 
In particular, the nucleon-nucleon interactions are normally not used directly 
in the nuclear many-body problem. Instead, effective interactions either "de- 
rived" from or motivated by the fundamental interactions are normally used 
in the many-body context f2lfl . Hence such interactions must necessarily be 
intimately linked to the Hilbert space of the calculation. Perhaps the most 
important ingredient of the many-body systems is to introduce the concept 
of the mean field. It is often sufficiently accurate to obtain a valid solution of 
the many-body problem. Based on the mean field assumption, the many body 
Schrodinger equation for the nucleus can be rewritten as 



6 



#Va(1,2,---,A)= ($> + ^ E %) V>a(1,2,---,A) (1) 
\i=i ^ ^=1 / 

= (H + >a(1, 2, • ■ ■ , A) = E^ x (l, 2, ■ ■ ■ , A), (2) 



with 



A 

HoMl, 2, • • • , A) = £ (ti + Ui ) 2, ■ • • , A) = EMI, 2, ■ • ■ , A) (3) 

%-X>), ( 4 ) 

where tj is the kinetic energy, and Ui is the mean field that the i-th nucleon feels 
in the many-body system, i.e., ui represents an average potential contributed 
by all nucleons in a many-nucleon system. H' is the residual interaction. 

Basically, there are two ways to deal with the above many-body equations: One 
is to focus on determining an optimal single-particle mean field Hq, namely, 
Hartree-Fock method. The other is to concentrate on figuring out a realistic 
two-body residual interaction H' in order to reproduce the spectroscopy of 
the many-body system. In the latter case, the technique assumes the solution 
to be a configuration mixture of Slater determinants. The solution divides 
neatly into two parts: the calculation of few body terms (one- and two- for 
most interactions) in a harmonic oscillator basis, followed by the many-body 
solution in second quantized form of eq.(|3]). This is the conventional Shell- 
Model approach. Prior to performing the Shell-Model calculation, one needs 
to derive the microscopic two-body interactions for a given model space. 

Both the Hartree-Fock and the Shell-Model approaches can also be used in 
conjunction with phenomenological interactions. In the first approach, the use 
of the Skyrme interaction leads to manageable calculations. It was shown that 



this interaction is related to the fundamental interaction ||22j| . On the other 
hand, for the Shell-Model, an approach based on a fit of the matrix elements in 
second quantized form @ attemps to differentiate the quality (or lack thereof) 
of the fit to experimental data as coming from either the many body aspects 
of the calculations or from the quality of the two body interaction. These 
approaches are fundamental in many aspects. 



4 Effective Nuclear Interaction for Shell-Model Calculations 



In the 60 's, Kuo and Brown presented a classic example to pertubatively de- 
rive the microscopic effective interaction (i.e., Shell-Model reaction matrix 
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Fig. 2. Comparison of microscopic two-body matrix elements in the s-d shell. 

elements) [p3|| . Later, a more exact and systematic way to derive the effec- 
tive two-body interaction for a given model space was developed, namely the 



Folded-diagram method. This method is reviewed in [24]. There is also a brief 
review on the NN effective interactions in [Ed* . 



For example, using the folded-diagram method, one can determine the 63 two- 
body matrix elements in the s-d shell, which together with the three single- 
particle energies are the essential input in the s-d shell microscopic calculation. 
Fig. 2 gives a comparison of those 63 matrix elements from Hamada- Johnston, 
Paris and Bonn-A potentials. It shows quite vividly that the different AW 
potentials can provide very consistent 2-body matrix elements. Only for some 
matrix elements, especially in the T = case, are there some discrapencies. 
However, these small discrapencies can produce substancial differences in the 
binding enengy calculations as we will see in the next section. 



5 Conventional Shell-Model Calculations 



Nuclei beyond the S-D shell require new techniques to handle the very large 
model spaces involved. Using massive amount of CPU on the fastest com- 
puters [pj| while keeping the traditional algorithms is simply not sufficient to 
model large nuclei. For this reason, there is renewed interest in recent years to 
develop better, faster and more robust Shell-Model algorithms. It deserves to 
be mentioned in this direction are the pair Shell- Model algorithm [27] and the 
Monte-Carlo approaches PB| . The former implements an exact Shell-Model 
solution for even-even nuclei within a truncated Shell-Model space built from 
arbitrary pair structures. It generalizes the FDUO code |29j code which was 
built from pairs envisioned in the Fermion Dynamical Symmetry Model 



30 



These codes can not address directly the question of the worthiness of the fun- 
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Fig. 3. Two coupling schemes for multishell calculations. 

damental interaction since they seek a solution in a severely truncated model 
space, which in itself would require an effective interaction with uncertainties. 
The Monte-Carlo Shell-Model code on the other hand seeks a solution of the 
many body problem via a Monte-Carlo variational approach. It appears to 
work will in reproducing the ground state properties of nuclei but suffers, at 
least up to now, from a lack of convergence for realistic interactions 



The Drexel University Shell-Model (DUSM) |31[] is another recently completed 
code to perform Shell-Model calculations which embodies an entirely new al- 
gorithm based on permutation group concepts. This is the code we used to 
perform the calculations to be described in this report despites their simplicity. 
This code promisses great performance in cases where a full exact Shell- Model 
solutions are seeked for even-even, even-odd and odd-odd nuclei. It outper- 
forms standard Shell-Model codes |$2] by large ( 4 or 5) factors in CPU usage 
with much reduced disk space and I/O requirements for J-T (spin-isospin) 
coupled spaces. 



5.1 Drexel University Shell Model Code (DUSM) 



The DUSM algorithm starts from the simple observation that Shell-Model 
calculations are often "multi-shell" in only a single subspace. For instance, a 
"J — T" calculation is multi-shell in the J subspace but single shell in the 
T subspace (t=| for all nucleons). Or a Shell- Model calculation for an elec- 
tronic system is multi-shell in L and single-shell in S. The wavefunctions in 
each subspace have arbitrary permutational symmetry; these symmetry pat- 
terns are simply restricted to be conjugate, so as to eventually couple to total 
antisymmetric wavefunctions. The DUSM coupling scheme is shown in the 
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grayed area of Fig. 3. The DUSM algorithm then proceeds to calculate wave- 
functions, coupling coefficents and matrix elements of elementary operators 
in each shell for each subspace and then in the total space according to the 
following steps: 



(i) Compute Coefficients of Fractioal Parentage (CFP) for each shell in each 
subspace by diagonalysing the Casimir operator of SU (2j + 1) in product 
basis. 

(ii) Compute the matrix elements in all single states of the elementary op- 
erators (all combinations of second quantized creation and annihilation 
operators) 

(iii) Compute the coupling coefficients for the coupling among shells (Outer 
Product Isoscalar Factors - OISF) 

(iv) Compute the coupling coefficients for the coupling between subspaces 
(Inner Product Isoscalar Factors - IISF) 

(v) Compute the matrix elements of the Hamiltonian via a "Sum over Path" 
in permutation diagrams concepts 

(vi) Diagonalyze the hamiltonian matrix 

(vii) Loop over all total quantum number (J and T for instance) to obtain the 
global solution 



The algorithm presents many advantages over traditional schemes: in partic- 
ular it avoids computing and storing total space CFPs. This implies much 
reduced I/O during execution and much less disk storage. Single shell calcu- 
lations are known to be possible up to very large J and number of particle; 
DUSM inherits this advantage. The DUSM algorithm uses group theory con- 
cepts at all levels to obtain the coupling coefficients, diagonalising the matrix 
represention of casimir operators of the appropriate groups at all steps, im- 
plying very stable numerical schemes. This is an advantage over using Racah 
formulas to derive the CFP [R3j in that no explicit orthogonalization is ever 



needed. The algorithm can furthermore be implemented (with pointers) with- 
out any search since the approach specifies fully the range over all intermediate 
sums. The approach is fully coupled; it provides full spectra and transitions 
among the states if required. 



This code has recently been completed; we use it in this report to compute 
some simple cases appropriate to understand the physics of the NN interac- 
tion. This does not illustrate by any means the capabilities of DUSM. 



5.2 s-d Shell-Model Calculations 



We now describe calculations done in various model spaces via DUSM using 
different effective interactions. 
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An important factor in the comparison of the experimental binding energies 
with results from large scale Shell-Model calculations lies in the evaluation 
of the Coulomb energy; this energy is subtracted from the total experimental 
binding energy. The relation between total binding energy and the Coulomb 
term Eq is as follow: 



Btotal — B 

core 

+ B SM + E* M , (5) 
E^ M = E l S tal - E c c ore , (6) 

where B core is the binding energy for the core. Bsm is the Shell-Model bind- 
ing energy related to the core; this term contains the binding energy of the 
cluster valence-particle plus the binding between the cluster and core. B c is 
the Coulomb energy. 

The total binding energy and core binding energy are listed in the 1993 atomic 
mass evaluation [Q. Therefore, according to eq. (f|), the binding energy cal- 



culated from Shell-Model can be directly compared with experimental values 
Bsm when the Coulomb term E^ M is extracted from the experimental values. 
In ref.[^3|, the Coulomb energy is estimated with two assumptions: first, for 
nuclei with a given proton number Z, the Coulomb energies are independent 
of mass, 

E C (Z,A) = E C (Z,A') . (7) 



Secondly, for nuclei with same mass A, the Shell-Model binding energies mea- 
sured with respect to the analogue isospin states are the same due to the 
isospin symmetry for the interactions between identical nucleons (i.e., proton- 
proton force and neutron- neutron force); consequently, Bsm for mirror nuclei 
are the same, while it differs from odd-odd nuclei due to the excitation energy 
measured with respect to the ground state, 

E J C T (Z,A) = E J C T (A-Z,A) . (8) 



For instance, when 16 is considered as an inert core, BsM{ 18 0)=BsM( 18 Ne) 
(the ground state is (JT)=01), while for 18 E, the excitation energy of the 
(JT)=01 state with respect to its ground state is 1.042 (MeV), thus B SM ( 1S F) 
+ 1.042 = Bsm{ 18 0). Notice that this way to estimate the Coulomb energies 
is not unique. Ref.[|35| claimed that the energy difference obtained via different 
routes is < 150/ceV under the above two assumption. 

In more realistic formulation, the Coulomb energy should be closely related 
to the proton distribution, which in turn should be affected by the neutron 
distribution. To be more specific, for example, with more neutron outside the 
core, the attractive interaction between proton and neutron tends to pull the 
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Table 2 

The empirical Shell-Model binding energies of ground states. 



(MeV) 


18q 


18 p 


18 Ne 


WB 


BsM 


-12.188 


-13.23 


-12.188 




E SM 





3.48 


7.666 


FRDM 1 


BsM 


-11.9855 


-12.7966 


-11.2471 




fiSM 


-0.2025 


3.0466 


6.7251 



protons towards the surface (i.e., disperse the charge distribution), thus reduc- 
ing the Coulomb energy. Therefore for nuclei with the same Z but different 
neutron number N, their Eq should be different. In fact, such effect is expected 
to be crucial for neutron rich nuclei. 



Alternatively, one can estimate the average Coulomb energy using an empirical 
formula |36| 



£ c = 0.717— (1 - —) (MeV) 



A 2 / 3 ' 



(9) 



The second term in equation above is the Coulomb exchange term which takes 
into account the effect of dispersion of nucleon due to Pauli exclusive effect. 
We mention that the theoretical results from the Fermi gas model agree nicely 
with this empirical formula. Recently, the finite-range droplet model (FRDM) 



| 36| was used to successfully evaluate nuclear ground-state masses and other 



nuclear structure properties^], including nuclei between the proton neutron 



drip lines, showing that the error did not increase with distance from (3 sta- 
bility. Similar to the Strutinski's method, FRDM describes nuclear binding 
energy with two components, one is the so-called macroscopic energy, which 
is,in some sense, taking care of nuclear bulk behavior, and the second com- 
ponent models the microscopic energy correction accounting for the nuclear 
shell structures and two-body residual interaction (mainly the pairing effect). 
The finite range effect of the nuclear force is accounted for by the macroscopic 
part. Thus the Coulomb energy is further modified [ 36| . In table 2, we give 



the empirical Shell-Model binding energies of 18 0, 18 F and 18 Ne using differ- 
ent methods p5| , |36| . Note that there are negative values of Eq , which means 



that Ec( 18 0) < Ec( 16 0) due to charge polarization. 

In Table 3, we list the single particle energies for A=18 nuclei. With the 
s.p. energies and 2-body matrix elements given from sec. 4, we now carry out 
the Shell-Model calculation for the binding energies of ground states. Table 4 
presents the calculated binding energies from the microscopic two-body matrix 
elements of Hamada- Johnston (i.e., Kuo-Brown (KB) interactions which were 
derived perturbatively in P3|), Bonn- A and Paris potentials. Comparing with 
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Table 3 

The s.p. energies of A=18 nuclei. 



(MeV) 


dr, 

2 


S 1 

2 


ds 

2 


18(9 E c =0 


-4.144 


-3.273 


0.941 


£ c =-0.103 


-4.041 


-3.170 


1.044 


18 F E c =3.543 


-4.144 


-3.649 


0.856 


E c =3.179 


-3.780 


-3.285 


1.220 




Fig. 4. Spectra of 18 0. 

Table 2, one can see that KB reproduces correct binding energies, However, 
Bonn results overbind about 1.8 MeV for T = 1 and 3.5 MeV for T = 0. Paris 
results are between KB and Bonn. Notice that the 2-body matrix elements 
from Bonn and Paris potentials are derived through a more precise method, 
namely the folded-diagram approach. These results of Shell- Model calculations 
are consistent with the few-body calculations with three-body force, where 
Bonn potential also gives rise to an overbinding. 

The above Shell-Model results for the binding energies are at best confusing: 
when one starts from a NN soft-core potential, based on effective meson the- 
ory with coupling constants being determined by analyzing about a few thou- 
sands of NN scattering data values and the deuteron properties, the binding 
energy results of the microscopic Shell-Model calculations are overbinding. 
On the other hand, starting from a NN hard-core potential, which is given 
perturbatively (i.e., the early KB effective interactions), the binding energy is 
reproduced very well. It is not only that the the binding energy is somewhat 
better, but the resulting spectra also demonstrate that KB's presents the best 
agreement compared with results from other modern potentials. We illustrate 
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Fig. 5. Spectra of 18 F. 



Table 4 

Comparison of binding energies 



Nuclei 




KB (MeV) 


Bonn (MeV) 


Paris (MeV) 


18 


kuo 


-12.20 


-13.96 


-13.02 




E c =0 


-12.19 


-13.95 


-13.01 




E c =-0.103 


-11.98 


-13.75 


-13.02 


18 F 


kuo 


-13.36 


-16.86 


-15.44 




£ c =3.543 


-13.58 


-17.09 


-15.68 




E C =3.179 


-12.85 


-16.36 


-14.95 


18 Ne 




-12.20 


-13.96 


-13.02 




E c =3.543 


-12.28 


-14.07 


-13.12 




E C =3.179 


-11.55 


-13.34 


-12.40 



this point in Fig.4 and Fig.5 where we show the theoretical and experimental 
spectra of 18 and 18 F. It is worth mentioning that KB not only provides the 
best microscopic shell model results for nuclei in s-d shell; recently, system- 
atic Shell-Model calculations for f-p nuclei have been carried out by using the 
Kuo-Brown interaction for the f-p shell PBJ showing a similar trend. 



We may point out that a main difference between the early KB effective inter- 
actions and the more recent ones (for example, see [p5|l ) is about the treatment 
of the folded diagrams. For the former the folded diagrams were ignored, with 
the effective interaction given merely by the bare-G and the second-order core 
polarization diagrams. These diagrams are usually referred to as G and G^ p ih 
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in the literature [37|. For the later, certain types of folded diagrams are in- 



cluded to all orders using a Q-box formulation |37j . What we have found in this 
work may indicate the need of a further investigation of the folded diagrams. 
There may be some other physical processes, which may counterbalance the 
effect of the folded diagrams and which have not been investigated. When this 
is done, it may be possible that the effective interaction will again be mainly 
given by G and G$ p \h alone, making life simpler. 



6 Comments 



In summary, we have given a brief review of the microscopic Shell-Model 
studies, and mentioned the recently developed Drexel University Shell Model 
(DUSM) code, which implements a new Shell-Model algorithm. We presented 
simple Shell-Model calculations to illustrate the progress made in this field 
and discussed some significant problems which remain in connection to the 
effective interactions used. 
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